Theory of orientational ordering in colloidal molecular crystals 
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Freezing of charged colloids on square or triangular two-dimensional periodic substrates has been 
recently shown to realize a rich variety of orientational orders. We propose a theoretical framework 
to analyze the corresponding structures. A fundamental ingredient is that a non spherical charged 
object in an electrolyte creates a screened electrostatic potential that is anisotropic at any distance. 
Our approach is in excellent agreement with the known experimental and numerical results, and 
explains in simple terms the reentrant orientational melting observed in these so called colloidal 
molecular crystals. We also investigate the case of a rectangular periodic substrate and predict an 
unusual phase transition between orientationnaly ordered states, as the aspect ratio of the unit cell 
is changed. 

PACS numbers: 82.70.Dd,64.70.Dv 



The theoretical scenario put forward by Kosterlitz, 
Thouless, Halperin, Nelson and Young for 2D melting 
has fostered many experimental investigations Most 
of these studies focused on 2D ordering on homogeneous 
substrates. On the other hand, the situation of a corru- 
gated substrate has received much less attention in spite 
of its experimental relevance, e.g to describe an atomic 
adsorbate on a crystalline surface. Colloidal particles 
subjected to light lattices created by laser arrays, appear 
as an ideal system to investigate the latter problem, see 
Refs I^LIlIIlI- These studies -mostly with ID periodic 
substrates- have revealed a novel laser induced freezing 
and laser induced melting that seems now well under- 
stood . It was only very recently that colloidal crystal- 
lization on 2D (triangular or square) periodic substrates 
has been investigated both experimentally Q and numer- 
ically @. The compatible results reported there call for 
a theoretical description that remains hitherto open. A 
rich variety of novel colloidal states were reported. When 
there are more colloids (number N c ) than the N substrate 
minima, the N c /N colloids trapped in each well may be 
considered as a single bound entity with only rotational 
degrees of freedom. These "colloidal molecules" (dimers 
for N c = 2N, trimers for N c = 3N. . . ) were found ex- 
perimentally and numerically to display long range orien- 
tational order. The corresponding state has been coined 
"colloidal molecular crystal" in Ref. Quite an un- 

usual two stage melting process may be observed for 
such systems. Decreasing the temperature T from an ini- 
tial standard 2D liquid, individual diffusion is first inhib- 
ited. The particles localize in the substrate wells (traps) 
thereby forming the colloidal molecules. The orienta- 
tional order in this partially ordered solid is short range. 
Upon further decreasing T, the colloidal molecules no 
longer rotate freely and acquire orientational order. Sur- 
prisingly, increasing the substrate strength Vq at fixed T 
leads to reentrant melting: the colloidal molecular crys- 
tal is destabilized into the partially ordered solid, and 
orientational order is lost. 

The aim of the present Letter is to provide a simple 
theoretical framework for the phase behavior of such a 



system, that is very different from that on homogeneous 
or ID periodic substrate potentials. We adopt here the 
simplest viewpoint and consider a trapped molecule (ri- 
mer) as a single rigid object with n-fold symmetry and 
one orientational degree of freedom within a common 
given plane. Additional degrees of freedom such as ra- 
dial or angular fluctuations will be discarded. We shall 
consequently not investigate the liquid phase (at high 
T and low Vq where the n-mers are dissolved) and con- 
centrate on the orientational freezing transition. In con- 
nection with numerical and experimental outputs, three 
questions naturally arise: a) is it possible to predict the 
non trivial ground state observed in the experiments and 
simulations? b) what are the characteristics of the above 
phase transition, in particular, its order? c) can we ex- 
plain the reentrant melting considering rotational degrees 
of freedom only? In addition, we provide predictions for 
further experimentally relevant situations. 



Before addressing these questions, general considera- 
tions concerning screened electrostatic interactions be- 
tween anisotropic charged molecules (the n-mers) are in 
order. Because their surface groups dissociate in solu- 
tion, the colloids under scrutiny here are highly charged. 
They strongly repel each other and experience a screened 
Coulomb (Yukawa) potential <j>(v) cx exp(— nr)/r, 
where the r dependence only is relevant for the discus- 
sion, not the prefactor depending both on particle size 
and charge |8( . In this expression, r is the center to cen- 
ter distance and n~ x the Debye screening length (smaller 
than the periodicity of the underlying 2D substrate). 
Once an n-mer is formed in a trap, how does it interact 
with a neighboring n-mer? Consider dimers for simplic- 
ity. One could naively think that in an electrolyte, the 
electrostatic potential created by such an object becomes 
isotropic at large distances, as is the case in vacuum or 
in a plain dielectric medium. This is however incorrect 
and summing the screened contributions from each parti- 
cle forming dimers i and j, one obtains the leading order 
large distance expression (see Fig. ^ f° r the notations 



2 



used) 

e~ Kr 

4>ij{r, 0ij,9ji) oc cosh[/tdcos(0y)] cosh[ndcos{9 ji)] 

r (l) 

The key point here is that the radial and angular de- 
pendence of the interaction potential are factorized. The 
anisotropy of a dimer is therefore felt at all distances 
even for «x ^> 1, provided the Debye length is small 
enough to resolve the dimer structure (when d <C l//t, 
the isotropic screened Coulomb potential is recovered). 
It is noteworthy to mention in passing that considering 
dipoles (two point charges +q/~q at a distance 2c?) in- 
stead of a dimer +q/+q, the same potential would hold 
provided the cosh function in Eq. is replaced by a 
sinh. The angular dependence of such an expression dif- 
fers much from that of the unscreened dipolar interaction. 
In addition, the distance dependences of both dimer and 
dipole potentials are identical (exp(— nr)/r in the far field 
region) . The same holds for higher order multipoles, that 
all contribute to the same order in distance. As a conse- 
quence, when considering the interactions between arbi- 
trary charged objects in a electrolyte, a usual truncated 
multipolar-like expansion |10| is doomed to fail; in the 
case of dimers, resuming all multipole orders would lead 
to Eq. to leading order in r. These considerations, 
that seem to have been overlooked in soft matter litera- 
ture, appear to account in simple terms for orientational 
freezing in colloidal molecular crystals. 




FIG. 1: Definition of the notations adopted to compute dimer 
(left) and trimer (right) interactions. The black dots show the 
colloids, and the line between them is a guide to the eye. 

In Ref. Reichhardt and Olson showed numerically 
that dimers on a square 2D periodic substrate (hereafter 
referred to as the underlying lattice of traps) adopt the 
ground state represented in Fig. 0-a). Such an order- 
ing should follow from minimizing the total electrostatic 
energy £ where each pair of dimers interact through the 
potential QJ. Given the exponential screening with dis- 
tance and that is smaller than lattice spacing, we 
further restrict to nearest neighbors of dimers and write 
the relevant angular dependent part of the energy as 

£ = cosh[Kdcos(Oij)) cosh[K<icos(#ji)], (2) 

(id) 

where the angular brackets denote nearest neighbors. A 
given pair minimizes its repulsion by setting 9ij = 9ji = 
ir/2 (parallel dimers, perpendicular to the center to cen- 
ter vector, see Fig. Such a pair configuration is how- 
ever not space-filling, so that no trivial ground state may 



be identified at this stage. By analogy with spin systems, 
we may invoke frustration. To proceed further, it is in- 
structive to approximate the cosh function in @ by a 
parabola. Introducing next a new set of variables where 
Oi G [— 7r/2,7r/2] is defined as the angle between dimer 
axis and one of the principal lattice directions ^2 we S e t 

£ = est + (kg!) 4 ( cose if ( cosd i) 2 ■ ( 3 ) 

(hi) 

Remarkably, this expression takes the form of an anti- 
ferromagnetic Ising-like Hamiltonian, once the variables 
<7i = cos26*i have been introduced (—1 < <Ji < 1): 

£ = est + {ndf ^criVj. (4) 

With the "spin" variables, frustration vanishes and the 
ground state simply consists of = ±1, with opposite 
signs on neighboring sites. This corresponds to 9i = 
or 7r/2, which is in perfect agreement with the results 
of the numerical simulations of Ref. [8j (see Fig. 0). 
We have checked that keeping the original expression @ 
without the parabolic approximation for the cosh leads 
to the same results. 



a) b) 




c) , , dl 



FIG. 2: Configuration of colloids (black dots) at T = on 
a rectangular lattice of traps, when iV c = 2N and nd — 1.5. 
a) situation on the square lattice (aspect ratio a = 1); b) 
rectangular lattice with a — 0.90; c) a = 0.85; d) a = 0.7 i.e. 
smaller than the threshold a* reported in Fig|3] Plots b) and 
c) correspond to the transition region a* < a < a** between 
F and AF phases. 

Having answered question a) above in the case of 
dimers on the square lattice, we now address ques- 
tions b) and c), in particular the reentrant melting phe- 
nomenon. At high temperature, the system has to be 
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in a disordered state, a paramagnetic phase in spin lan- 
guage, and on purely dimensional grounds, Eq indi- 
cates that the critical temperature of the antiferromag- 
netic/paramagnetic transition scales like T c oc (nd) 4 . At 
fixed temperature T, increasing the substrate pinning 
strength Vq, the n-mers constituting a colloidal molecule 
become of smaller spatial extent: d decreases, so that 
T c eventually becomes smaller than T. Orientational or- 
der is correspondingly lost. More specifically, we find 
T c cx V Q 4 ^ 3 . As emphasized above, we discard the possi- 
ble loss of positional order, that occurs at large T. The 
corresponding fusion temperature Tf has to increase with 
Vq. For large enough Vn, the transition between the 
ordered and partially ordered solids is therefore always 
present, since T c < Tf. At low Vo, this transition is pre- 
empted by the positional melting (Tf < T c ). 

The above arguments may be rationalized by comput- 
ing the partition function associated to the Hamiltonian 
(I2J) UM- A mean-field analysis, simplified by the intro- 
duction of two sub-lattices to describe the order displayed 
in Fig. leads to the conclusion that orientational freez- 
ing is a second order phase transition |14| . 

The same procedure may be extended to other sub- 
strate geometries, or more complex colloidal molecules. 
For n-mers, the counterpart of Eq. Q reads 



Ojj(r.O, ,.Q ,,) = f(9ij,Kd)f(9ji,Kd)' 



where 



n-l 

f(8, Kd) — exp 

k=0 



ndcos 8 



2kir\ 
n J 



(5) 



(6) 



2). 



which reduces to Eq. JQ) in the case of dimers (n 
The variable 9ij is defined as the angle between one arm 
of the symmetric rt-fold colloidal molecule i and the cen- 
ter to center direction (ij), see Fig. \\\ior the cases n = 2 
and n = 3. For the situation of a triangular substrate 
with trimers as colloidal molecules (the experimental sit- 
uation investigated in |7]]), we may follow similar steps 
as above, i.e. minimize ^ f(@ij)f(@ji) where the sum- 
mation is restricted to nearest neighbors. This leads to 
the ground state where all trimers point in the same 
direction, with an angle 8 = 7r/3 between one arm of 
the trimers and one of the principal lattice direction (or 
equivalently, fly = 7r/3 ± 7r/3 depending on the pair (i- 
j)). This is exactly the order observed experimentally 
Q, and in numerical simulations For dimers on 

a triangular substrate (not studied experimentally), we 
also obtain the same orientational order as in Ref. Q 
(herringbone structure with 8 — ±7r/4, or equivalently 
fly =7r/4±7r/3). 

Our framework allows to investigate the situation of 
non-integer filling. We will focus on the half integer case 
N c /N = 3/2, as considered numerically in ^{|. To mini- 
mize electrostatic repulsion on the square lattice, the col- 
loids create a checkerboard order with two sub-lattices 



such that the nearest neighbors of monomer are exclu- 
sively dimers, and vice- versa [l5| . In this configuration, 
nearest neighbors interaction do not depend on the an- 
gular coordinates, so that the electrostatic energy to con- 
sider should include second nearest neighbor terms. We 
find that the corresponding ground state is degenerate, 
with dimer orientations 8i — or 8i = 7r/2. Since there 
are N c dimers, this ground state is 2 Nc degenerate, and 
no long range order can be observed. This finding sup- 
ports the numerical results reported in [l5| . Changing 
the geometry of the underlying substrate from square to 
triangular does not allow to tile the lattice such that min- 
ima with only one colloid will be in between every other 
dimer . The partitioning of the 3N C colloids in the 27V 
traps is therefore a frustrated problem with no ordering, 
be it positional or orientational. 




FIG. 3: Phase diagram for dimers on a rectangular substrate 
as a function of the aspect ratio a. The F phase is obtained for 
a < a* , while a > a** corresponds to the AF region. Here 
nd is fixed to 1.5, which would correspond experimentally 
to a fixed depth Vo- For a given value kI — 6 (horizontal 
dashed line in the main graph such that a* ~ 0.82 and a** ~ 
0.93), the inset shows how the system evolves from the F 
to AF phase as a increases. The cross, triangle and square 
correspond respectively to the situations of Fig. |5Jd) with 
a = 0.7, Fig. 0c) with a = 0.85, and Fig. 0b) with a = 0.90. 



We now consider a final geometry of experimental rele- 
vance, namely dimers on a rectangular lattice (with unit 
cell of lengths al and I; by definition, a < 1, and the di- 
rection along the long axis is called "horizontal" for con- 
venience) . When the aspect ratio a is 1 the checkerboard 
order is that of Fig. 0-a) (antiferromagnetic-like phase, 
denoted AF). On the other hand, when a<l, there are 
only two nearest neighbors per particle, and from expres- 
sion JQ| the preferred configuration is fly = Bji = ir/2: 
all dimers are parallel to the horizontal direction. This 
ferromagnetic-like phase is denoted F. We have studied 
the transition between the F and AF phases, minimizing 
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the relevant energy 

£rcct = Pij cosh.[Kdcos(8 i j)]cosh[Kdcos(8j i )), (7) 

where /3y = a exp[/tZ(l — a)] (resp. /3y = 1) for a 
vertical (resp. horizontal) bond i-j. As expected, when 
a is close to 1 (in fact a > a**), the order is the AF 
phase of Fig. 0-a) and when a is small enough (a < a*), 
the F phase of Fig. EJ-d) is observed. As shown in Fig. 
13 the thresholds a* and a** depend on screening condi- 
tions and dimer length through nl and nd. Investigating 
the behavior in the transition region a* < a < a** re- 
quires the inclusion of next nearest neighbors |l4j . and 
leads to the conclusion that the system divides into two 
sub-lattices A and B (see Fig|2-b) and c)). On A, the 
dimers remain parallel to the horizontal axis while on 
B, they have alternating orientations ±# depending on 
the row (see Fig. |5J). The inset of Fig shows how the 
orientations evolve from the F state with 9 = to the 
AF order with 9 = tt/2. The F/AF transition driven by 
a is therefore second order. We have also studied more 
complex situations with higher fillings (integer and half 
integer), that realize a host of non trivial colloidal crys- 
talline states flij . 

Our approach relies on the validity of the pair poten- 
tial Q which is an expansion obtained at large distances 
from the original Yukawa expression. Such an expansion 
is justified for d <C I and (nd) 2 <C kI. In addition, we 
restricted to nearest or second nearest neighbors interac- 



tions, which requires kI ~^$> 1. Knowledge of the precise 
parameter bounds where our predictions apply -in par- 
ticular on the triangular lattice- is an important question 
of experimental relevance, that can only be answered nu- 
merically. Equally relevant is the influence of tempera- 
ture, that may affect the scenarios reported here for the 
ground state behaviour (such as those of Figs and |2J . 
Work along these lines is in progress, with Monte Carlo 
simulations. 

In conclusion, our theoretical approach explains in sim- 
ple terms the rich orientational phenomenology of col- 
loidal molecular crystals, both the emergence of long 
range order and reentrance. Our work establishes a 
natural connection between the situation of dimers on 
a square lattice and the canonical Ising model of sta- 
tistical mechanics. The possibility of such a mapping 
however seems lost for more complex colloidal molecules 
or with other substrate geometries. The mechanisms at 
work here discard the radial and angular fluctuations of 
a trapped colloidal molecule, that appear to be of sec- 
ondary importance. On the other hand, such fluctuations 
are prevalent on one-dimensional periodic substrates [(| 
and could a priori have played a role in the present two- 
dimensional situation Q- Finally, we hope that our re- 
sults will stimulate further experimental investigations in 
this field. 

We would like to thank C. Bechinger, S. Bleil, L. Boc- 
quet, M. Brunner, B. Jancovici and Z. Racz for useful 
discussions. 
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